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Abstract 

The fc-space polarization structure and its strain response in SrTi03 with rotational instability 
are studied using a combination of first-principles density functional calculations, modern theory 
of polarization, and analytic Wannier-function formulation. (1) As one outcome of this study, 
we rigorously prove — both numerically and analytically — that folding effect exists in polarization 
structure. (2) After eliminating the folding effect, we find that the polarization structure for SrTi03 
with rotational instability is still considerably different from that for non-rotational SrTi03, re- 
vealing that polarization structure is sensitive to structure distortion of oxygen-octahedra rotation 
and promises to be an effective tool for studying material properties. (3) Furthermore, from po- 
larization structure we determine the microscopic Wannier-function interactions in SrTiOs. These 
interactions are found to vary significantly with and without oxygen-octahedra rotation. 

PACS numbers: 77.80.-e, 77.84.-s, 77.22.Ej 
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I. INTRODUCTION 



Materials with perovskite or perovskite-based structures have attracted wide attention 
for their varieties of properties in ferroelectricity, magnetism, superconductivity, etc. In 
these perovskite solids, an important class of structural distortions that deviate them from 
ideal (cubic) perovskites is the rotation of the oxygen octahedra, which may considerably 
alter the physical and chemical properties of materials. For instance, rotation of oxygen 
octahedra in incipient ferroelectric SrTiOs was found to couple strongly with tetragonal 
strain, and consequently, affect significantly the structural and dielectric behaviors below 
105K. Octahedral rotation also plays a key role in improper ferroelectrics such as 

YMnC>3, where structural instability due to zone-boundary phonon mode couples with elec- 
tric polarization [3], and this coupling was recently shown to exhibit an interesting absence 
of critical thickness in improper-ferroelectric ultrathin films. |4| On magnetism, Mazin and 
Singh have revealed that one main reason responsible for the drastically different magnetic 
properties in two chemically similar compounds CaRu0 3 and SrRu0 3 (where CaRu0 3 is 
paramagnetic, but SrRu0 3 is ferromagnetic) consists in a larger rotation angle of RuOg 
octahedra in CaRu0 3 .|l]] Similarly, Fang and Terakura found that perovskite superconduc- 
tor Sr2Ru04 could be turned from non-magnetic to ferromagnetic by simply increasing the 
rotation angle of the RuO@. f| Recently, Ray and Waghmare have shown that rotation is 
the predominant structural instability in a series of biferroics (LuCrOa, YCrOq, LaCrOs, 
BiCrOs, and YCrOa), and is closely correlated with their magnetic properties. [7| Based on 
the facts that (1) rotation of oxygen octahedra affects various material properties, and (2) 
the oxygen-octahedra rotation is common, understanding of the physics caused by oxygen 
rotation is thus important. 

The modern theory of polarizationjs, 9] has become important in quantitatively determin- 
ing the electric and dielectric properties of materials within first-principles density functional 
calculations. According to this theoryjsl, 0], the change of electric polarization is a physi- 
cally meaningful quantity and the electronic contribution is determined by an integration 
of the geometric phase of the occupied multi-fold electron Bloch wavefunctions. Based on 
the modern theory of polarization, Yao and Fu recently introduced the concept of "polar- 



ization structure" 



10[ which may provide a new scheme to study the electronic properties of 



dielectric materials, and they went on to formulate a theory using localized Wannier func- 
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tions to investigate the polarization dispersion structure of ferroelectrics with spontaneous 
polarization. More specifically, polarization structure describes the electronic phase angle <fi 
as a function of k± point which lies on a /c-plane perpendicular to the direction of the con- 
sidered polarization component. At each k±, the phase angle is determined by the electron 
wavefunctions asjs, 9] 

M f G \\ „ d 



71=1 



where subscripts || and _L indicate, respectively, the directions parallel and perpendicular 
to the direction of the polarization. u n ^ is the Bloch part of electron wave function, and M 
is the number of occupied bands. As significant as electronic band structure, polarization 
structure reveals how individual strings of different fcj_s contribute to the electronic polar- 
ization, which thus offers microscopic insight on the polarization properties. Polarization 
structure can be changed by applied electric fields, in a similar spirit as electronic state in 
band structure can be changed by photoexcitation. More discussions on the importance 



of the polarization structure were given in RefJlOl. It was also shown that the polarization 
structure is determined by the interaction among Wannier functions in neighboring cells. [h3] 
Polarization structure of a dielectric system may thus enable scientists to study the micro- 
scopic interaction among Wannier functions in real space, which is of fundamental relevance 
in revealing material properties. 

In this work we address a topic which concerns both the rotational instability of oxygen 
octahedra and polarization structure, by investigating how the polarization dispersion struc- 
ture in ferroelectric perovskites is influenced by oxygen rotation. We will use SrTi03 (ST) 
as an example. ST is known as an incipient ferroelectric, which exhibits no polarization 
at the strain-free ground state but is ready to develop a polarization under small external 
strains. UJ It was observed experimentally that unstrained ground state of ST has a nonzero 



oxygen-octahedra rotation around the (001) direction. Later, we will see from our numerical 
simulations that by applying in-plane compressive strain, the rotation angle in ST increases 
and the structure with oxygen rotation remains to be the most stable ground state as the 
strain is small or moderate. But at large strains, the tetragonal phase of zero rotation be- 
comes more stable. Comparing SrTiOa (with oxygen rotation) with PbTiOs (of no oxygen 
rotation), one anticipates both the atom-atom interaction and electron wave functions to be 
altered in two materials. Since the polarization structure is determined by electron wave- 
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functions and their interaction, it is thus of interest to study how the polarization structure 
changes as a result of the octahedral rotation. It worths pointing out that the influence of 
rotational instability on the polarization structure has not been studied before 10] , and turns 
out to offer interesting new knowledge on the phase contribution at individual kj_ points as 
well as on microscopic interaction. 

The purpose of this work is three folds: (i) To determine and reveal how the polarization 
structure may look like in ferroelectrics with rotation instability, using first-principles density 
functional calculations; (ii) By comparing the polarization structure in SrTiC>3 with oxygen 
rotation (labeled as R-SrTi0 3 ) with the counterpart in non-rotational SrTi0 3 (to be denoted 
as NR-SrTi03 where the oxygen-octahedral rotation is deliberately suppressed by symmetry 
constraint), we intend to investigate the effect of rotational instability on the formation of 
polarization dispersion structure; (iii) Since octahedral rotation couples with, and can be 
changed by, inplane strains, we thus would like to explore how the polarization structure in 
perovskites with rotation instability may be modified by the inplane strain. 

The content of this paper is organized as following. In Sec. II we describe the computa- 
tional methods used in this study, along with a brief result on the structural responses in 
R-SrTiC>3 and NR-SrTiCV In Sec. Ill we present the polarization structures for R-SrTiC>3 
and NR-SrTiC>3, and identify the differences and similarities between these polarization 
structures. Two main reasons contributing to the differences are revealed. In addition, 
analytical formulation is also made in Sec. Ill to deduce from polarization dispersion the 
microscopic interaction among Wannier functions. We find that the 3rd nearest neighbor 
(3NN) approximation is necessary in order to yield an accurate polarization structure for 
NR-SrTiC>3 under large inplane strains. Based on our analytical formulation and numerical 
simulation, we also describe a renormalization-group scheme to determine the higher-order 
interactions between Wannier functions. Sec. IV summarizes the results of this work. 



II. METHODS AND STRUCTURAL RESPONSE 

We use first-principles density-functional theory (DFT) within local density approxima- 
tion (LDA) [12[ to determine cell parameters and internal atomic positions in zero-strained 
SrTiC>3 as well as in SrTiC>3 under different compressive in-plane strains. More specifically, 
for a given inplane lattice constant a, both the out-of-plane lattice constant c and atomic 
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locations are optimized by minimizing the total energy and forces. For the optimization of 
atomic positions, residual forces on each atom are required to be smaller than 10 -3 eV/A. 
Effects of nuclei and core electrons are replaced by Troullier-Martins norm-conserving type 
of pseudopotentials [13 ] ; semi-core states Ti 3s and 3p are treated as valence states for better 
accuracy. Bloch wave functions in solid are expanded using a mixed-basis set consisting of 
numerical atomic pseudo-orbitals and plane waves. 14j A cutoff energy of 100 Ryd in wave- 
function expansion is used for both rotational and non-rotational SrTi03 on equal footing, 
and has been checked to be sufficient. Details on generating pseudopotential, pseudo/all- 



electron matching radii, and the accuracy checking were described in Ref. 



15) . Polarization 



structure, namely the contribution of a given kj_ string to the electronic polarization, is cal- 



culated according to the equation 



M 0(^0 = Im{ln n j=0 det((u mifcj > nifcj+1 ))}, where J is 
the number of k\\ points along the string. One advantage of this equation over Eq.((T]) is that 
(p(k±) here does not depend on the arbitrariness of the phases of Bloch u n ^ wavefunctions 
at individual k point. 

Figure [Ha) and (b) show the top- view of the cell structure for NR-SrTi03 without octa- 
hedral rotation of 5-atoms each cell, and for R-SrTi03 with octahedral rotation of 10-atoms 
each cell, respectively. In NR-SrTi03, the lattice vectors are a'^ai, a' 2 =aj, a' 3 =ck. Biaxial 
compressive strains are applied in the a\ and a! 2 directions. In R-SrTi03, the lattice vectors 
are di=a{i — j), a 2 =a(z + j), a 3 =ck. The rotation angle 9, as illustrated in Figfjjb), char- 
acterizes the deviation of the oxygen atoms from the non-rotational positions in the a\-a2 
plane. In the rest of this section, we will briefly give our results on the optimized struc- 
tures for NR-SrTi03 and R-SrTi03 under different inplane strains, since these data could 
be useful for experiments or other theoretical calculations. The main results on polarization 
dispersion will be discussed in the next section. 

Figure |2]^a) gives the total energies per 10 atoms for NR-SrTiC>3 and R-SrTiC>3 at different 
in-plane a lattice constants; each curve is obtained by imposing symmetry constraint on the 
corresponding structure. One main result in Figj2]Ja) is that a crossover in the ground-state 
structure occurs at a c =3.62A. When a > a c , R-SrTi03 has a lower energy than NR-SrTi03 
and is thus more stable. As a is below a c , NR-SrTi03 becomes more stable, and the crystal 
of strained SrTi03 changes from rotational to non-rotational. Meanwhile, we need to point 
out that at the equilibrium lattice constant a =3.86A, our LDA simulations yield a rather 
small energy difference between NR-SrTi03 and R-SrTi03, which is on the order of numerical 



5 



uncertainty, and does not allow us to tell which structure is more stable. This, however, does 
not contradict with the fact that R-SrTi03 is found more stable in experiments, since LDA 
calculations can not capture the zero-point quantum fluctuation which is known important in 
stabilizing the experimentally-observed rotational ground state of incipient SrTi0 3 . On the 
other hand, since the main purpose of this work is to understand the polarization-structure 
difference between R-SrTi03 and NR-SrTi03, the above uncertainty in the ground-state 
structure does not affect our purpose of investigating the polarization structure. 

The calculated c/a tetragonality is given in Table HJ along with the rotation angle 9 in 
R-SrTi0 3 . Under zero strain, both R-SrTi0 3 and NR-SrTi0 3 are found to possess an equi- 
librium inplane lattice constant of ao=3-86A, agreeing well with the experimental value 16] of 
3.90A for R-SrTi0 3 at low temperature. For R-SrTi0 3 at a , our theoretical rotation angle is 
#=3.25°, comparable to the experimental value |l7j of 2.1° as well as 5.5° in another theoret- 



ical calculation 



181 ]. Under inplane strains, the c/a ratios are shown in FigfS^b). One could 



see that by decreasing a, c/a increases for both NR-SrTi0 3 and R-SrTi0 3 . At small strains 
(a>3.7lA), c/a's of two structures remain close; however, at large strains (a<3.7lA), c/a 
of NR-SrTi0 3 increases dramatically, and becomes evidently larger than that in R-SrTi0 3 . 
Interestingly, for the strain range studied here, c/a of R-SrTi0 3 depends on a in a rather 
linear fashion. In contrast, for NR-SrTi0 3 , the c/a ratio undergoes an obvious slope change 
around a=3.7lA. Strain-induced total polarization is shown in Figj2]^c), where we see (i) for 
NR-SrTi0 3 , there is a critical inplane lattice constant a m 3.82A, above which polarization 
remains null; (ii) However, for R-SrTi0 3 , such a critical inplane lattice constant apparently 
does not exist, and instead polarization increases almost linearly with the decreasing inplane 
lattice constant. 

The differences between NR-SrTi0 3 and R-SrTi0 3 in the c/a tetragonality [Figj2](b)] and 
in the electric polarization [Figj2]^c)] can be explained microscopically by the ferroelectric 
atomic off-center displacements, which are shown in FigG^d). In NR-SrTi0 3 , atomic dis- 
placements of 01 and 02 atoms are small when a>3.80A and then increases drastically 
with a notable slope change as a is below 3.80A, which causes the sharp rise of electric 
polarization in Figj2]^c). On the other hand, in R-SrTi0 3 the off-center displacements of 01 
and 02 vary rather linearly with the decreasing inplane lattice constant, without significant 
change in slope. At a=3.59A, the oxygen displacements in NR-SrTi0 3 are more than two 
times as those in R-SrTi0 3 . 
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III. POLARIZATION STRUCTURE 



A. Comparison on polarization dispersions 

We now present and discuss the polarization structures for SrTiC>3 with and without 
oxygen-octahedra rotation, in particular the similarity and difference between them. Under 
zero strain, both NR-SrTi03 and R-SrTi03 have null polarizations, and the 4>{k±) phase 
is trivial and vanishes for all k± points. On the other hand, under compressive strains, 
both systems develop electric polarizations. In the following, we chose an inplane lattice 
constant a=3.71A, and determine the polarization structure for NR-SrTi03 and R-SrTiC>3 
by using the optimized atomic positions and cell shapes of each system. Since NR-SrTi0 3 
and R-SrTi03 have different sizes of unit cells, their Brillouin zones (BZ) are illustrated in 
FiglS^a) and (b), respectively. The triangles highlighted by thick lines in Figj3]^a) and (b) are 
the high-symmetry paths along which the /cj_-dependent polarization structure is examined. 
The calculated <f)(k±) phase structures are shown in Figj3]^c). To facilitate comparison, both 
polarization structures are shifted to make 4>(T) = at T; meanwhile, the horizontal axis 
of k± distance in R-SrTi0 3 is rescaled by a factor of \/2 in Figj3]^c), such that the two 
polarization dispersions can be given in the same figure. 

Inspection of the polarization structures in FigfS^c) reveals: (i) The dispersion bandwidth 
in NR-SrTi0 3 (~0.9) is considerably larger than that in R-SrTi0 3 (~0.4). (ii) In NR-SrTi0 3 , 
the largest (f>(k±) contribution comes from the X2 point at the zone corner, while in R-SrTiC>3 
the largest contribution is at Mi which is on the zone boundary, (iii) In NR-SrTiC>3, Xi 
point is a saddle point, while in R-SrTiC>3 there is no saddle point. Furthermore, in R- 
SrTiC>3, M 2 is a local minimum and another local maximum appears along the M 2 — T line, 
which is not the case in NR-SrTiC>3. As a result, the two polarization dispersion curves have 
rather different curvatures. 

While the above differences in the polarization structure between NR-SrTiC>3 and R- 
SrTi03 are interesting, they are also puzzling. With some effort, we find two possible reasons 
that might be responsible for these differences. The first one is folding effect. It is known in 
band-structure calculations that when one changes the size of unit cell, there is band-folding 
effect. However, it is unclear whether the folding effect also exists for polarization dispersion. 
If the answer is yes, one may also desire to establish a rigorous proof on the existence of 
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such an effect, and better yet, to uncover the precise manner on how the 4>{k±) dispersion 
curve folds itself. The second possible reason is rotational effect. When octahedral rotation 
is allowed, the atom-atom interaction in R-SrTi03 is anticipated to be different from in 
NR-SrTi03, which could alter the electron distribution and thus the polarization structure. 
In the following, we will determine and analyze the consequences of both effects. 



B. Folding effect 

To determine the significance of each of the above two effects (i.e., folding effect and 
rotational effect), it is preferable to design an approach that is able to investigate each effect 
separately. We first investigate the folding effect, by conducting the following calculations 
in order to avoid deliberately the rotational effect so that the latter will not play a role. 
We begin with an LDA-optimized 5-atom cell for NR-SrTi03 at a given inplane lattice (the 
solid of which is to be denoted as NR5), and then double the size of cell and convert it into 
a 10-atom cell (the solid of which is to be denoted as NR10, once again without oxygen 
rotation). The atomic positions and the electron charge densities in NR5 before conversion 
and in NR10 after conversion are thus equivalent. Also note that there is no oxygen rotation 
in either NR5 or NR10; differences in polarization structure between these two solids thus 
come from the folding effect alone. 

FigJH(a) presents the polarization structures of NR5 and NR10, both at a=3.7lA. One sees 
that the polarization structure of NR10 (empty triangles in FigH^) appears very different 
from that of NR5 (empty squares). In particular, (i) The bandwidth of the NR5 curve 
is about twice as large as the NR10 case, (ii) For NR5, X2 is the maximum, and X\ is a 
saddle point as we saw previously. However for NR10, M\ and M2 have similar 4>(k±) values, 
making the curve between Mi and M2 evidently flat, (iii) The curvatures are different in 
two solids. As a result, apparently no relationship could be found for the two polarization- 
dispersion curves at first glance. More disturbingly, one may also wonder whether the 
polarization structure is uniquely defined for a given crystal, or whether it is meaningful to 
make a sensible comparison on polarization dispersion. To answer this, we will establish a 
rigorous relationship between the two curves in FigJUa), using an analytical approach based 
on Wannier functions. 



— * 

We start with the expression for the polarization phase structure [lOj, cj){k 
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T" J2r ± J2n=i J r\\W*(r)W n (r — R±)e tkj -' R± df, where W n (r) is the Wannier function of the 
solid defined as W n {f — R) = J BZ dke' lk ^ r ~ R ^u n k(r) ■ Subscript || (or JL) in the equation 
means the projection component parallel (or perpendicular, respectively) to the direction of 
the considered polarization; R± thus refers to a lattice vector perpendicular to the direction 
of the polarization. M is the number of occupied bands in solid. By defining parameters 

t(R±) = — E / nWn(?)W n (r-R ± )dr, (2) 

° n=l J 

we then arrive at a concise and analytic expression for polarization structure, 

^k ± ) = J2t(R±)e rh ^ . (3) 



Hi 



— * 

Sum over R± in Eq. ([3]) should in principle been carried out to infinite distance. In 



19 



practice, 



20|, it is 



owing to the localization nature of Wannier functions in dielectric materials 
often sufficient to truncate R± within a few near neighbors (NNs). 

To simplify the formulations, we will consider, for NR5, up to the 3rd nearest-neighbor 
interactions in Eq.Q, with parameters t , t 1; t 2 , and t 3 representing respectively the on-site, 
first NN (INN), second NN (2NN), and third NN (3NN) Wannier-function interactions as 
schematically illustrated in FigfjTJc). On an equal footing, the same truncation distance 
is used for Wannier-function interactions in NR10 where cell size is doubled. However, 
we should point out that our final conclusion on folding effect is generally applicable and 
does not depend on the distance within which Wannier-function interactions are truncated. 
Furthermore, to avoid ambiguity, we will use R ; ± and <j)' for NR5, while we use R± and 
for NR10. For NR5 within 3NN truncation, R' ± s are R' ± =(0,0) (on site); (±a,0) and (0,±a) 
(INNs); (±a,±a) (2NNs); (±2a,0) and (0,±2a) (3NNs). With tetragonal symmetry, we can 
rewrite Eq.Q as 

0'(^_l) — + 2t 1 [cos(fcia) + cos(fc 2 a)] 

(4) 

+4t 2 cos(fcia) cos^a) + 2t 3 [cos(2/c 1 a) + cos(2A;2a)] , 

where k\ and k 2 are the components of kj_, namely k± = (k\, k 2 )- 

For NR10, there are two sets of Wannier functions in one unit cell, one set being the 
same as those in NR5 and another set having identical shapes but displaced by the lattice 
vector a[ of NR5. Meanwhile, we denote the Wannier-interaction parameters in NR10 as Tj 
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as schematically shown in FigJT^d), to be distinguished from tj in NR5. Within the same 
truncation distance as in NR5, we only need to consider in NR10 the interactions between 
Wannier functions up to 2NNs, where R± = (0,0) (on site); (±a,±a) (INNs); (±2a,0) and 
(0,±2a) (2NNs). The polarization structure in NR10 could be determined as 

<t>{k±) = T + 4Ti cos(fcia) cos(A; 2 a) 

(5) 

+2T 2 [cos(2fci<z) + cos(2fc 2 a)] . 
By the definition of Eq.([2]), it can be further shown that £jS in NR5 and TjS in NR10 are 
related by To = 2to, T\ = 2t 2 , and T 2 = 2t%, since the Wannier functions are identical in 
both systems. Combining Eq.(j4]) with Eq.flSJ), one can find the following relation between </> 
and (j)' phases 

<t>(k±) = <f/(k±) + <f>'(k±-b 2 ) , (6) 

where b 2 = -{i + j) is one of the reciprocal lattice vectors of NR10. In other words, for an 
arbitrary k± point within the triangle TMiM 2 in the irreducible BZ of NR10, the phase angle 
4>{k±) of NR10 could be determined from the phases of NR5, at two k± points within the 
triangle rXxX 2 in the irreducible BZ — one contribution being at the same k± point (the first 
term in Eqj6]) and the other is the folded contribution at the k± — b 2 point (the second term 
in Eqj6]). We thus prove that the folding effect indeed exists for the polarization dispersion 
structure. 

We next would like to numerically confirm the validity of Eq.((6]). We apply Eq.([6]) and 
use the right-hand side of this equation to determine the polarization structure <f)(k±) of 
NR10. The results are given as solid squares in FigJD^a), where we see that the analytic 
results of Eq.()6]) coincide with those LDA calculated results of NR10. This demonstrates 
that the polarization structures of NR5 and NR10, although quite different looking, are 
equivalent and are thus physically meaningful. 

C. Rotational Effect 

After revealing the existence of folding effect, we proceed to investigate rotational effect. 
To single out the rotation effect, we compare in Figj4(b) the polarization structures of 
NR10 (without rotation) with that of R-SrTiC>3 (with rotation), both at the inplane lattice 
constant a=3.71A. Note that both systems have 10 atoms in one unit cell, and the folding 
effect is irrelevant here. 
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After eliminating the folding effect, we see in FigH^b) that the bandwidths of the po- 
larization structures with and without rotation become notably similar. However, the two 
polarization structures still show important differences: (i) 0(M 2 ) is a global maximum in 
non-rotational NR10, whereas being a local minimum in rotational structure; (ii) <f>(Mi) is 
still a saddle point in NR10 as in NR5, but is a global maximum in rotational case; (iii) 
Between Mi and M 2 , very little dispersion exists for the polarization structure of NR10, 
which is not the case for R-SrTi03 where sizeable dispersion is evident; (iv) A local maxi- 
mum is observed along M 2 — T in R-SrTi03, and such a maximum does not exist in the case 
of NR10. Rotation of the oxygen octahedra alters the interactions among atoms (also the 
Wannier functions), which is responsible for the above- described differences in polarization 
structure. Differences (i)-(iv) also reveal that polarization structure is sensitive to the subtle 
changes in atomic positions and crystal structure, and could be used as an effective tool for 
studying material properties. 

D. Response of polarization structure to in-plane strain 

Compressive inplane strains — caused by lattice mismatch between epitaxially grown solid 
and substrate — modify interatomic interactions in an anisotropic fashion, by strengthening 
the interaction along inplane directions but weakening the interaction along the out-of-plane 
direction. This modification affects the structural and electric properties. As a consequence 
of the altered interatomic interaction, the polarization structure might also change accord- 
ingly. In this part, we examine the strain dependence of the polarization structure in R- 
SrTiOa with octahedral rotation, and contrast it with that in the non-rotational structure. 
To eliminate the folding effect, we use a unit cell of 10 atoms for non-rotational SrTiOs (i.e., 
NR10), same as for R-SrTiC>3. The k± points are along the T — > Mi — > M 2 — » T path as in 
FigE^b). The calculated polarization structures at different in-plane lattice constants are 
shown in Figj5l More specifically, for R-SrTiC>3, its polarization dispersions under different 
inplane a lattice constants are shown in Figj5]^a), while Figj5]^b) summarizes the phase 4>(k±) 
angles at high-symmetric Mi and M 2 points as a function of the inplane lattice constant. 
For non-rotational NR10 structure, the polarization dispersions are given in FigH^c) while 
phase angles at Mi and M 2 are summarized in FigJS^d). 

Let us focus on R-SrTiC>3 first. In Figl5](a), we see that with the increase of in-plane 
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strain, the polarization structures change as following: (i) At zero strain or a=3. 86A, the 
polarization <f)(kj_) phase vanishes at all k± points, as a result of the existing inversion 
symmetry; (ii) Upon application of small strain, the bandwidth of polarization structure 
begins to increase. For 3.83A > a > 3.71A, the 4>(k±) angle peaks at M 1; while being 
a local minimum at M 2 . Meanwhile, a local maximum is observed along M 2 -T. On the 
other hand, the zone center, T, contributes the least to the total electronic polarization; (iii) 
When 3.7lA > a > 3.65A, 0(M 2 ) increases while </>(Mi) slightly decreases, making the (j)(k±) 
curve of a=3.65A notably flat between M\ and M 2 . (iv) As a is further decreased below 
3.65A, M 2 turns into the maximum point, while Mi becomes a saddle point. Meanwhile, 
the local maximum previously existing along M 2 -T disappears. Quantitatively, evolutions 
of the 4>{k±) phases at Mi and M 2 with the inplane lattice constant are better visualized 
in Figj5]^b). The evolution can be separated into three regions: Region I where <fr(Mi) and 
<f)(M 2 ) both increase with the increasing strain, Region II where 4>(Mi) starts to decline but 
4>(Mi) still contributes greater than <f)(M 2 ), and Region III where crossover between <j)(Mi) 
and 0(M 2 ) occurs and <fi{M 2 ) now becomes the maximum contribution. 

The differences between R-SrTiC>3 and NR10, in their responses of polarization structure 
to inplane strain, can be seen by comparing FigJS^d) with Figl^b), which yields the following 
observations: (i) When a is near 3.83A, 0(Mi) and 0(M 2 ) both increase linearly in R- 
SrTiC>3 [Figj5](b)], while they remain nearly zero in NR10 [Figj5]^d)]. It reveals once again a 
lack of critical inplane strain in R-SrTiC>3. However, in NR10 without rotation, a need be 
below 3.83A in order to produce sizable polarization, (ii) When 3.83A< a < 3.74A, 4>(M 2 ) 
increases quickly with strain both in R-SrTiC>3 and in NR10, but the increasing slope in 
NR10 [FigJS](d)] is much larger than in R-SrTiC>3 [Figl^b)]. During this region of strain, 
the polarization in NR10 is more sensitive to the inplane strain with larger piezoelectric 
response than in R-SrTiC>3, and microscopically <fi{M 2 ) is largely responsible, (iii) Although 
the responses in Figjf^d) can also be divided into three stages I, II and III for NR10, the 
inplane lattice constants delimiting the boundary of different stages differ, however, from 
those in FigJS](b). (iv) In NR10, 4>(Mi) declines sharply at stages II and III, and as a result, 
4>{M 2 ) becomes dominating when a < 3.68A. In R-SrTiC>3, 0(M X ) declines much slowly at 
stages II and III. (v) 0(M 2 ) in NR10 reaches its maximum around a=3.65A, then starts to 
decrease. In R-SrTiC>3, <p(M 2 ) nevertheless keeps increasing. 

Response of polarization structure to inplane strain is thus found interestingly different 
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in R-SrTi03 with oxygen rotation as compared to that in NR10 without rotation. These dif- 
ferences should, on the microscopic level, originate from the different interaction among the 
Wannier functions in two structures. In the following, we will extract from the polarization 
structure the interactions among Wannier functions. 

E. Implication on interaction between Wannier functions 

Our procedure of determining Wannier-function interaction begins with Eq.flS]), which was 
derived for crystals consisting of 10 atoms in one unit cell, and is thus applicable to both R- 
SrTiC>3 and NR10 though with different interaction parameters. Applying Eq.Q to <p(k_i_) at 
special k± points at T, Mi and M 2 , we obtain </>(Mi)—</>(T) = — 4Ti— 8T 2 and 0(M 2 ) — (f>(T) = 
—8T\. From <p{Mi) and 0(M 2 ), one can thus determine 7\ and T 2 quantities. As shown in 
Figfjjd), T\ (or T 2 ) corresponds to the Wannier-function interaction between a given cell with 
the 1st (or 2nd) neighboring cell, manifesting interatomic interaction on microscopic level. 
The obtained 7\ and T 2 quantities are then used in Eq.(j5]) to determine analytically the whole 
polarization dispersion curve. Note that only phases at M\ and M 2 points are fitted. To 
evaluate the validity of our analytic formulas, we choose three inplane lattice constants across 
the entire strain range considered in this work, and plot the analytic polarization dispersions 
as open symbols in Figl^a) for R-SrTiC>3 and in FigJH^b) for NR10. One sees that for each 
inplane lattice constant, the polarization structures analytically determined from Eq.([I)]), and 
the direct results obtained from the Berry-phase modern theory of polarization, agree well. 
More specifically, for the a=3.80A curve in Figlf^a), the small dispersion hump between 
M 2 and V is well reproduced by analytical result. Furthermore, at a=3.59A which is a 
more stringent test on the analytic formulation with truncated interaction, the saddle-point 
nature at Mi in FigJI^a), as well as the subtle minimum at Mi in FigJS](b), are also predicted 
by the analytic expressions. These demonstrate the validity of our analytic formulation; the 
obtained T\ and T 2 values are thus reliable. 

The obtained T\ and T 2 quantities are given in Table HT1 and shown in Figj6^c) and (d). 
One can see from Figlf^c) and (d) that, for a given inplane lattice constant, Tj is notably 
different in R-SrTi03 with respect to that in non-rotational NR10. For instance, at a=3.68A, 
Ti in NR-SrTi03 is ~50% larger than in R-SrTi03, and T 2 in NR-SrTi03 even has a different 
sign from T 2 in R-SrTi03. These quantitative differences in microscopic interactions indicate 
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that polarization structure is indeed a rather powerful tool to understand physical properties. 
We also note that, at large inplane strain such as a=3.59A, T 2 is large in NR-SrTi03 but 
nearly vanishes in R-SrTi03 [FigJBfd)], implying that interactions in rotational and non- 
rotational structures are significantly different. Furthermore, we observe that at very small 
inplane lattice constants near a=3.59A, T\ and T 2 tend to saturate in NR-SrTi03 but not 
in R-SrTi03. 

F. Determination of higher-order Wannier-function interaction 

One relevant question in the theory of polarization structure is the determination of 
the tiS quantities, particularly those higher-order £jS that correspond to larger R±s, since 
by knowing t^s the <p(k±) phases can be readily obtained via Eq.(j3]). The £, quantities 
connect the microscopic interaction with macroscopic fcj_-dependent polarization structure. 
As one useful outcome of the current study, we describe an approach that can be utilized to 
determine the higher-order t( interactions. 

The approach works in general. But for the sake of convenience in our discussion, 
we will take NR-SrTiC>3 (without oxygen rotation) of a=3.59A as an example. For NR- 
SrTiC>3, the smallest unit cell consists of 5 atoms, and the polarization structure directly 
calculated from Berry phase is given in FigfTl Now if we use a lower-order 2NN ap- 
proximation in Eq. (f3"]) , with interactions to, t\ and t 2 in Figfjjc), to analytically mimic 
the polarization structure, it is easy to derive that the polarization structure will be 
4>(k±) = to + 2ti[cos(fcia) + cos(fc 2 a)] + 2t 2 [cos(fci + k 2 )a + cos(A;i — ^2)0]. By using the 
phase angles at T, X\ and X 2 points, one can determine to, ti and t 2 . The polarization 
dispersion within the 2NN approximation can be subsequently obtained and is given in 
FigJTl One sees that the agreement between the directly calculated polarization curve and 
the 2NN-approximated curve is not good. For instance, 4>{X 2 ) is incorrectly predicted by 
the 2NN approximation to be a local minimum, instead of a local maximum as by the direct 
calculation. Meanwhile, the dispersion curvatures along X\-X 2 and along X 2 — V are wrong. 
We thus see that higher-order interactions among Wannier functions are needed. This is not 
surprising, since at small inplane lattice constants, interactions between farther neighboring 
cells remain strong and can not be neglected. 

To determine higher-order tj quantities, one possible approach (denoted as approach I) 
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is to use higher-order expression such as Eq.fll]) within the 3NN approximation, and 

in Eq.pj) four tiS need be fitted. In this approach, if even higher-order approximations than 
3NN are used, more t{S have to be fitted at many k± points, which becomes increasingly 
complex and less trackable. An alternative approach (denoted as approach II) is the fol- 
lowing. First, we recognize from the finding of the folding effect in Sec.IIIB that, when 
one doubles the size of unit cell from 5-atoms to 10- atoms, the solid is the same, but the 
t\ term drops as seen by comparing Eq.Q with Eq.(j2]). As a result, one is still left with 
three (not four) parameters, and can be determined by fitting 0(fc_i_)s at three k± points. 
Similarly, one can determine even higher-order Wannier-function interactions by varying 
the size of unit cell. In fact, when one increases the size of unit cell (doubling, tripling, 
etc), the interactions between those Wannier functions — that are located in different cells 
of previously smaller cell size but fall in the same cell of the larger cell size — are dropped. 
As a result, with the large unit cell, determination of new Wannier-function interactions 
becomes simplified, which corresponds to higher-order NN interactions with respect to the 
previous smaller unit cell. For instance, in Eq.([5]) when cell size is doubled, a new interaction 
parameter T 2 (i.e., £3 since T 2 = 2^) is determined. This approach is, in spirit, similar to 
the renormalization-group (RG) approach used in understanding critical behaviors of phase 



transition 



2l| . Of course, the benefit of approach II is at the expense of computing the 



Berry-phase polarization structure in larger unit cell, but the increase of computation is 
rather small. 

To examine how approach II works, we give in Fig J7] the analytical dispersion determined 
within the 3NN approximation, with t\ obtained from 5-atom-cell calculation, and with £ 2 
and £3 obtained from 10-atom-cell calculation. The obtained values are ti=-3.590xl0~ 2 , 
t 2 =-2.711x 10~ 2 , and t 3 =1.761x 10 -2 . As one could easily see from FigJTJ inclusion of t 3 
interaction in the 3NN approximation yields a much better dispersion than the 2NN one. The 
local- maximum nature at X 2 , the characteristic dip between X\ and X 2 , and the curvature 
between X2 and T in the direct Berry-phase calculation results are well described by the 
3NN approximation but not the 2NN approximation. 
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IV. SUMMARY 



Polarization structure in R-SrTi03 with rotational instability is studied, and contrasted 
with the dispersion in NR-SrTi03 where oxygen-octahedra rotation is deliberately not al- 
lowed by symmetry constraint. Influences of in-plane compressive strains on the polarization 
dispersion as well as on the phase contributions at high-symmetric k± points are investi- 
gated. Microscopic interactions are deduced by means of Wannier-f unction formulation, and 
the dependencies of T\ and T 2 interactions on the inplane strain are determined. 

Our specific findings are summarized as following: (i) The polarization structures of ST 
with and without oxygen rotation, at the same in-plane lattice constant, are rather different. 
These differences could be attributed to two reasons, the folding effect and the rotational 
effect, (ii) We have proved, both numerically and analytically via the Wannier-function 
formulation, that folding effect exists in polarization structure. As a result, polarization 
structures are equivalent for differently chosen unit cells, and are thus physically meaningful. 
We also reveal how polarization structure folds itself, (iii) After elimination of folding effect 
by choosing the same unit cell, the polarization structures for rotational and non-rotational 
geometries have similar bandwidth. But the dispersion curvatures and individual <f){kj_) 
phases are still considerably different as shown in FigJD^b), which reveals that polarization 
structure is sensitive to the structure distortion and is indeed an effective tool to study 
materials properties of dielectrics, (iv) With application of in-plane compressive strain, 
the <f){k±) phase at high-symmetric Mi point is found to decline sharply at stage III in 
NR-SrTi03 [FigJS^d)], and in contrast, the decline is much slower in R-SrTiC>3 [FigJS^c)]. 
Furthermore, in NR-SrTiC>3, the phase at another high-symmetric M 2 point first 

increases at small strain and then declines at high strain. However, in R-SrTiC>3, the <fi 
phase at M 2 keeps increasing throughout the considered strain range, (v) From polarization 
structure, we obtained the microscopic 7\ and T 2 interactions as given in Figj6]Jc) and (d). 
Ti interaction is found to differ significantly in R-SrTiC>3 and NR-SrTiC>3. Moreover, T 2 in 
Figlf^d) is revealed to exhibit an interesting non-monotonous dependence on inplane lattice 
constant, (vi) We have further described a renormalization-like method to determine higher- 
order Wannier-function interactions by varying the size of unit cell. We revealed that the 
3NN approximation is necessary for an accurate description of the polarization structure of 
NR-SrTi0 3 at a=3.59A. 
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FIG. 1: (Color online) (a) Top view along the [001] direction of the 5-atom cell, for SrTiO"3 without 
oxygen-octahedra rotation, (b) Top view of the 10-atom cell for SrTiO"3 with oxygen rotation, (c) 
Schematic illustration of the interaction between Wannier functions in different cells, for NR-SrTi03 
consisting of 5-atom cells, (d) Similar as (c), but for R-SrTi03 consisting of 10-atom cells. 
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FIG. 2: (Color Online) The following calculated properties as a function of the in-plane a lattice 
constant in R-SrTi03 and NR-SrTiOs: (a) total energy of a 10-atoms cell (with respect to a 
basis of -7474 eV), (b) the c/a tetragonality, (c) the total (electronic and ionic) polarization, (d) 
atomic displacements along the c-axis. Atomic displacements in (d) are obtained by comparing the 
LDA-optimized atomic positions with respect to the high-symmetric locations. High-symmetric 
locations are defined as: we first shift the Sr atom to the origin of the coordinate, and then the 
high-symmetric atomic positions are z = for Sr and 01 atoms, and z = | for Ti and 02 atoms. 
01 is apical while 02 is on the base plane of an octahedra. 
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FIG. 3: (a) The two-dimensional Brillouin zone over which k± is sampled, for NR-SrTi03 with 
five atoms in one cell, (b) The two-dimensional Brillouin zone over which k± is sampled, for R- 
SrTiOa with ten atoms in one cell, (c) Polarization structures of NR-SrTi03 and R-SrTi03, both 
at a=3.7lA. For NR-SrTi03, 5-atom cell is used for the polarization structure calculation. For 
R-SrTiO"3, the horizontal kj_ distance is rescaled by a factor of y/2 such that its dispersion structure 
and the counterpart of NR-SrTiO"3 can be plotted in the same figure. 
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FIG. 4: (a) Polarization structures of NR5 and NR10 at inplane lattice constant a=3.7lA. The 
analytic result using Eq.© is also shown as solid squares, (b) Comparison between the polarization 
structure of NR10 and that of R-SrTi03 at a=3.7lA. 
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FIG. 5: (Color online) (a) Polarization structure of R-SrTiO"3 at different in-plane lattice constants; 
(b) Phase angles at Mi and M2 as a function of the in-plane a lattice constant, for R-SrTiO"3. (c) 
Polarization structure of non-rotational NR10 at different in-plane lattice constants; (d) Phase 
angles at M\ and M2 as a function of the in-plane a lattice constant, for NR10. All polarization 
dispersion curves are rigidly shifted such that (j>(T) = for the sake of comparison. 
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FIG. 6: (a) Analytic polarization dispersions (open symbols) determined from Eq.([5]), as compared 
to the direct results (filled symbols) calculated from the modern theory of polarization, for R- 
SrTiO-3 at inplane lattice constants a=3.80, 3.65, and 3.59A. Lines are guides for eyes. Since 
analytic results agree well with direct calculation results, most of open symbols are on top of filled 
symbols, (b) Similar as (a), but for NR10 at a=3.77, 3.71, and 3.59A. (c) Dependence of the T\ 
quantity on the inplane a lattice constant, for R-SrTi03 and NR10. (d) Dependence of the T2 
quantity on the inplane lattice constant. 
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FIG. 7: Polarization structure of NR-SrTi03 at a=3.59A computed from direct Berry Phase calcu- 
lation (solid squares), 2NN approximation (open circles and dotted line), and 3NN approximation 
(open triangles). Lines are guides for eyes. 



TABLE I: The LDA-optimized c/a ratios for NR-SrTiC>3 and R-SrTiO"3 under different inplane 
lattice constants. For R-SrTiO"3, the calculated rotation angle is also given. 





NR-SrTi0 3 


R-SrTi0 3 


a(A) 


c/a 


c/a 


angle 0(°) 


3.59 


1.32 


1.21 


12.70 


3.62 


1.28 


1.17 


12.31 


3.65 


1.21 


1.14 


11.71 


3.68 


1.16 


1.12 


10.92 


3.71 


1.10 


1.09 


9.95 


3.74 


1.068 


1.067 


8.80 


3.77 


1.040 


1.050 


7.81 


3.80 


1.026 


1.031 


6.57 


3.83 


1.014 


1.016 


5.39 


3.86 


0.998 


1.002 


3.25 
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TABLE II: T\ and T 2 quantities (in radian) at different inplane lattice constants. 





NR10 




R-SrTi0 3 




a (k) 


T x (10~ 2 ) T 2 (10~ 2 ) 


Ti (10~ 2 ) T 2 (10~ 2 ) 


3 BQ 


-b.422 


o coo 

3.522 


coin 

-6.210 


-0.372 


3 62 


-o. i\)o 


Q CIO 

6.012 


C OC£ 


-1. / 06 


3.65 


-6 408 


9 849 


-4 848 


-9 714 


3.68 


-6.234 


1.350 


-4.064 


-3.116 


3.71 


-5.162 


-1.888 


-3.334 


-3.636 


3.74 


-4.058 


-3.612 


-2.780 


-3.662 


3.77 


-2.272 


-3.298 


-2.322 


-3.412 


3.80 


-0.476 


-0.778 


-1.552 


-2.562 


3.83 


9.74E-2 


1.74E-1 


-0.562 


-0.998 
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